Dense granular flows: 
two-particle argument accounts for friction-like constitutive law with threshold 
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A scalar constitutive law is obtained for dense granular fiows, both in the inertial regime where 
the grain inertia dominates, and in the viscous regime. Considering a pair of grains rather than a 
single grain, the classical arguments yield a constitutive law that exhibits a flow threshold expressed 
as a finite effective friction at flow onset. The value of the threshold is not predicted. The resulting 
law seems to be compatible with existing data, provided the saturation at high velocity (coUisional 
regime) is added empirically. The law is not exactly the same in both regimes, which seems to 
indicate that there is no "universal" law. 
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I. INTRODUCTION 

Dense flows of dry granular materials and granular 
pastes is still a surprising field of research Many 
complex features have already been observed, but the 
basic mecanisms that relate grain properties to the col- 
lective behavior are not fully identified yet. Gaining new 
insights into this multi-scale issue should provide robust 
models and rationalize a wide range of results obtained 
in various geometries and with various materials. 

Here, we focus on the central question of the consti- 
tutive law, which relates the stresses (shear stress r and 
pressure P) to the shear rate 7 in a steady and homo- 
geneous flow. Several measurements evidenced a flow 
threshold with a frictional criterion: a flow stops if the 
ratio between the shear stress and the pressure (effective 
friction coefficient ^ = -p) is lower than a critical friction 
coefficient fic- This is reminicent of yield stress fluids, 
which stop below a critical shear stress Tc. But for gran- 
ular materials this yield stress depends on the pressure 
(tc = fJ-cP). Measurements also evidenced an increase 
of the friction coefficient with the shear rate. Further- 
more, dimensional analysis implies that the effective fric- 
tion coefficient cannot depend just on the shear rate 7, 
but rather on a dimensionless combination / = Tj where 
T is some time scale. T is usually understood as the typ- 
ical time for a single grain, accelerated by the pressure 
P (or force P(P), to move over a distance comparable 
to its own size d. It has be estimated in two situations 
through a simple dimensional analysis. For dry grains, 
the surrouding fluid can be neglected and the pressure is 
resisted only by grain inertia (mass m) 0]: 

m^^Pd' ^.e., T,, = ^ — . (1) 
For granular pastes, if the grain inertia is negligible as 
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compared to the viscous effect of the surrounding fluid, 
the force Pd^ exerced by the pressure is resisted by a 
typical Stokes viscous drag force felt by a single grain 
moving at the velocity d/Tf^ [^: 

^d^^Pd^ i.e., Tfv = ^. (2) 

Using these time scales provided robust scaling behav- 
iors in various flow geomerties for dry grains, whether 
frictional or not 0, [j, H, @| , cohesive [3, 111 or polydis- 
persed Q, and also for immersed grains Empirical 
functional form for the constitutive law were proposed to 
describe the measurements 0, i, i, i, H, 0, S, H [IS [HI, 
[l3|,for instance (3]: 

This shape includes a critical friction coefficient below 
which flow stops, fic- At large /, a saturation of the 
effective friction (/x(/) — *■ ^2) reflects the onset of the 
collisional regime, where grains are diluted and interact 
through binary collisions [1, 0, H]. In the dense regime, 
at small /, the effective friction is assumed to increase 
linearly with / (slope ), and thus with the shear 

rate. However, a quadratic term in ^ was also included 
to fit the data in other works [13, El, So far, there 
is no consensus on the shape of the constitutive law. As 
for the choice of a frictional form like ([3]), it is justified 
purely phenomenologically. 

In this paper, we propose a simple analysis to derive 
such a frictional constitutive law from the grain proper- 
ties. We apply this analysis first to dry grains with no 
fluid, then to immersed grains with no inertia, and fl- 
nally to a situation with both the grain inertia and the 
fluid viscosity. We compare the model thus obtained to 
the experimental data with immersed and dry grains by 
Cassar et al. Q, kindly provided by the authors. 

Hopefully, the argument presented here may later 
be useful to understand the origin of much complex 
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FIG. 1: Schematic evolution within the granular material 
during shear, (a) One grain is transported quasistatically 
from position 1 to position 2, then falls into position 3 due 
to the applied pressure P. In the present work, we rather 
consider a pair of grains during the period of time when they 
are close neighbours (bl-3). First, the deviatoric (typically, 
shear) component of the stress, r, helps the pressure P es- 
tablish the contact (bl). Then, the contact rotates due to 
the overall material deformation (b2). Finally, the deviatoric 
stress overcomes the pressure to break the contact (bS). Be- 
cause the pressure is compressive in a non-cohesive granular 
material, the typical magnitude of the force transmitted be- 
tween both grains is stronger when the contact forms than 
when it breaks. 



features of dense granular flow: the viscosity bifurca- 
tions [H, [3, [i^ , a critical shear rate below which no ho- 
mogeneous flow exists 0, H, 0, [13, 13 ^ possible coexis- 
tence of flowing and non-flowin g region s (shear-banding, 
shear localization, cracks) 0, [tI ll7l.Tl8l. [TqI . [20| . 



II. ONE-GRAIN VERSUS TWO-GRAm 
ARGUMENT 

As mentioned above, the classical approach consists in 
assuming that the constitutive law has a frictional form 
like -p = and to include a characteristic time such 
as those provided by Eqs. (P) and that reflect the 
motion of a single grain subjected to the pressure P in 
the granular material. 

As shown on Fig. [l] the present approach is focused 
on the entire lifetime of a contact between two grains 
rather than on the motion of a single grain. The contact 
is beeing established when the pair of grains is oriented 
along a compressive direction, sec drawing (bl). Hence, 
the deviatoric part of the stress is complemented by the 
pressure part to favour the approach of both grains. The 
opposite is true when the contact is breaking (b3): the 
pair is now oriented along a tensile direction, and the 
pressure partly resists the deviatoric part of the stress. 
The normal force acting on the contact is thus stronger 



during the approach than during the separation: 

'd^P + r) 
(P - r) 
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(4) 
(5) 

(6) 

where r > is the magnitude of the deviatoric stress. 
Note that the stress on the contact evolves as described 
above even if the macroscopic stress does not vary. In- 
deed it is the rearrangement of the neighbouring grains 
that induce the progressive reorientation of the pair of 
interest and the corresponding inversion of projection of 
the deviatoric stress on the pair axis. 

Let us now assume that during a 100% deformation, 
all inter-grain contacts are typically renewed. The cor- 
responding time 1 /j is then comparable with the typical 
lifetime of a contact: 
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where the value of T^pp and T'^p is given by Eqs. dH) 
and replacing the pressure by r ± P as in Eqs. (|3][S]): 
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From Eqs. ([7]) and ^S^, wc obtain a local scalar con- 
stitutive equation for the inertial regime: 
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which can be approximated as: 
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since the first factor in Eq. varies smoothly between 
1 and 1/2 for t > P. 

Similarly, using Eqs. ([7]) and pOlllip . we obtain a local 
scalar constitutive equation for the viscous regime: 

^ ^'^ (14) 
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III. DISCUSSION AND COMPARISON WITH 
EXPERIMENTS 

When combined with Eq. (IT|) or ([2) as appropriate, 
Eqs. ([HI) and ([HI) become: 
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Even though these two equations do not reduce to 
Eq. ([3]), let us emphasize that they imply the existence 
of a flow threshold: flow can occur (7 > 0) only when 

T>P. 

However, the value of the threshold, named hi 
Eq. ([3]), is not predicted in the present approach. Indeed, 
the qualitative argument presented earlier, that leads to a 
different value for the force between two particles during 
their approach and during their separation is not precise 
enough to provide any reliable prefactor for r in Eqs. ([4]- 

ED. 

As shown by Fig. [21 the saturation of function 
at large / introduced empirically in Eq. ([3]) is not repro- 
duced by Eqs. (|15m6p . as the arguments used to obtain 
the adimensional parameter / for the dense flows are not 
compatible with the binary collisions thought to be at the 
origin of the saturation. However, the plots of Eqs. (1151 
I16p arc not very different from that of 
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(see curve Coo of Fig. [2]), obtained from Eq. Q with 
/ic = 1, after taking the limit fj.2 ^ 00 and Iq 00 while 
keeping Iq = ^2/Mc~1i in order to remove the saturation. 

Let us introduce empirically the very same saturation 
into Eqs. (|15m6p in order to be able to test whether 
they are compatible with the data presented by Cassar 
et al. 

Eq. ^ can be rewritten as: 
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Similarly, let us transform Eqs. (|15m6p into: 
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Eqs. p8ll20p are plotted on Fig. [3] in the form ^ as 
a function of /, together with the data of Fig. 13 from 
Rcf. 3, kindly provided by the authors. 

It appears that with adjusted values of /q and /i2, 
Eqs. (|19ll20p are also compatible with the data. More 
elaborate comparisons should be carried out to see 
whether the aerial and submarine data can be usefully 
distinguished using the above equations. 

For completeness, an interpolation between Eqs. (fTS]) 
and pop has been derived in the same manner and is 
available in the form of a short note [211 . 



IV. CONCLUSION 

Considering a pair of grains rather than a single grain 
in a dense granular flow suggests that the typical stress 




FIG. 2: Comparison of the present predictions with the em- 
pirical function /i(/). The curves labeled "gi" and "fv" are 
the plots of Eqs. (|15m6|l respectively. Curve "C" corresponds 



to Eq. (131, with /^c = 1, ^2 = 2 and 7o = 1. Curve 
obtained after removing the saturation, see Eq. (fT7)l . 
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FIG. 3: Effective friction coefficient /i(/): data by Cassar et 
al. 01 kindly provided by the authors, and plot of three differ- 
ent functions. Open circles (Wl): 112/^m diameter beads in 
water. Filled circles (Al): same beads in air. Open squares 
(W2): 208^m beads in water. Filled squares (A2): same 
beads in air. Curve "C" corresponds to Eq. ^ with Jo = 1, 
/ic = 0.43 and ^2 ~ 0.82. Curve "gi" represents Eq. (|19|l with 
/o = 0.2, /ic = 0.43 and /i2 — 0.7. Curve "fv" corresponds to 
Eq. ((201) with Jo = 0.7, /ic = 0.43 and /t2 = 1. 



during the approach of two grains differs from that during 
their seperation. This difference implies that the usual 
arguments used to derive an intrinsic timescale, whether 
when the grain intertia is dominant or when the fluid vis- 
cosity is dominant, now additionally predict the existence 
of a frictional threshold, with flow only when the ratio of 
the deviatoric stress to the pressure exceeds some finite 
value: t/P > Hc- However, the threshold value /ic (usu- 
ally around 0.4) cannot be predicted by this approach. 

Cassar et al. had obtained a "universal" behaviour 
in the empirical form of Eq. in terms of a dimcn- 
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sionless parameter / with a different meaning in botfi Acknowledgements 
regimes. Tlie fact tliat two distinct laws are found for 
the inertial regime and for the viscous regime with the 
present approach indicates that the universal behaviour 

is perhaps only approximate. More involved investiga- This work was supported by the Agence Nationale de 
tions will be needed to test this universality. la Recherche (ANR05). 
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